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Abstract 

The effect of quenched disorder in the one-dimensional asymmetric exclu- 
sion process is reviewed. Both particlewise and sitewise disorder generically 
induces phase separation in a range of densities. In the particlewise case the 
existence of stationary product measures in the homogeneous phase implies 
that the critical density can be computed exactly, while for sitewise disor- 
der only bounds are available. The coarsening of phase-separated domains 
starting from a homogeneous initial condition is addressed using scaling ar- 
guments and extremal statistics considerations. Some of these results have 
been obtained previously in the context of directed polymers subject to 
columnar disorder. 

1 Introduction 

The one-dimensional asymmetric simple exclusion process (ASEP) was introduced 
by Spitzer in 1970 as an example of an interacting stochastic process @, |<§. 
In the probabilistic community it has been widely used for rigorous studies of 
the emergence of hydrodynamic behavior from stochastic microscopic dynamics 



P, Already thirty years ago similar models were considered in the context of 
biopolymerization ||, while recent applications have focused on the problem of 



vehicular traffic flow [^]. The interest of statistical physicists has been further 
fueled by the discovery of boundary- induced phase transitions || |, [H| as well as 



the relations to interface growth and directed polymers in random media [11], [12 
In short, the ASEP is a generic model of driven single file transport which combines 
utmost simplicity with a remarkable richness of behaviors. 

Figure 1 illustrates the model. Particles occupy the sites of a one- dimensional 
lattice subject to the simple exclusion rule (at most one particle per site). In an 
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infinitesimal time interval dt particle % at site Xi attempts a jump to the right (left) 
with probability pdt (qdt). The jump succeeds if the neighboring site is empty and 
is suppressed otherwise. In general the jump rates p and q may depend on both 
the particle label i and the position x on the lattice. In much of the paper I will 
restrict myself to the totally asymmetric case q = 0. 

In the present article I want to address the effects that quenched disorder in 
the jump rates has on the behavior of the ASEP. Disorder effects can be quite 
dramatic in one-dimensional single file systems, as is evidenced by the everyday 
experience with platoons and traffic jams caused by slow vehicles, accidents or 
road construction on highways |13| , [14], 15], |16, 17[] . Also in the context of driven 



transport on biomolecules a certain amount of disorder seems unavoidable 

It is natural to distinguish between particlewise disorder with p = Pi-, q = q% 
independent of x, and sitewise disorder with p = p(x), q = q(x) independent of 
i. Both particlewise and sitewise disorder generically induces phase separation in 
the sense that, for global particle densities p in a certain interval [p~,p+], the 
system breaks up into regions of density p~ and p+ separated by sharp density 
discontinuities ("shocks"). These shocks are typically associated with bottlenecks, 
i.e. slow particles or slow sites in the particlewise and sitewise cases, respectively. 
If the system is started from a homogeneous initial condition, the average size £ 
of phase separated regions grows as a power law 

i{t) ~ fV', (1) 

defining a dynamic exponent z\ an example of the time evolution in the particlewise 
case is shown in Figure 2. Two kinds of questions will therefore be asked in 
the following: First, how can the density interval [p~ , p+] of phase separation be 
determined? Second, what is the value of the dynamic exponent, and how does it 
depend on the distribution of the disordered jump rates? 



For particlewise disorder a number of exact analytic results are available [|19 



p0| , |2T| , [22] , p3| which have been reviewed elsewhere |Tj|. This case will therefore 
only be briefly summarized in Section [|. The more difficult problem of sitewise 



disorder has been studied numerically by Tripathy and Barma [G3 and others 



jig , [25] , but little is known analytically. In Section y some progress in this di- 
rection will be reported. Specifically, I derive a rigorous bound on the critical 
densities based on the results for the particlewise case, and obtain predictions for 
the coarsening behavior for various types of disorder distributions. The relation to 
directed polymers in random media is briefly discussed in Section |3.4| , and some 
conclusions and open questions are formulated in Section [|. 



2 



2 Particlewise disorder 



2.1 Steady state and critical density 

For particlewise disorder the configurations of the system are most naturally de- 
scribed in terms of the headways = Xi + % — X{ — 1 in front of the particles. The 
key simplifying feature is that different headways become statistically independent 
in the steady state, with a geometric distribution |19|, |21[ 



Pi(u) = (1-eiiK (2) 

for the headway in front of particle i. In the totally asymmetric case the parameters 
ctj are determined by the jump rates pi through the simple relation p0| , ^H] 



on = v/pi (3) 

where v is the (common) mean speed of the particles in the steady state. Eq. @ 
expresses the plausible fact that the headways in front of slow particles are larger 
than in front of fast ones. The geometric distribution (|2|) remains valid in the 
case of partial asymmetry, but then (Q) is replaced by a more complicated relation 
|19| , pTf . The steady state distribution for the totally asymmetric model with 
parallel update has a similar form |22] . 



In the following we consider the totally asymmetric case and take the pi to be 
independent random variables with a probability density f(p) supported on the 
interval [c, 1], with a minimal speed c bounded away from zero. Since particles 
cannot pass each other, it is clear that the steady state speed v in an infinite 
system cannot exceed c. To compute it, one determines the mean headway in front 
of particle i from (0) and performs the disorder average. In a system of density p 
the resulting average headway must be (1 — p)/ p. This yields the implicit equation 



P 



1 + v 



1 dp f(p) 



P 



(4) 



for the speed as a function of density. Two cases are to be distinguished. If the 
integral on the right hand side of (Q) diverges in the limit v — > c, then v (p) < c 
for all p > 0. In this case the «j in @ are bounded away from unity for all i, the 
headway distributions are normalizable, and the system remains homogeneous. If, 
on the other hand, the integral remains finite in this limit, then the right hand 
side of (£|) evaluated at v = c defines a critical density p c such that v (p) = c 
in the entire interval [0,p c ]. For the slowest particles with pi w c this implies 
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that the headway distributions (|2|) are no longer normalizable. Large gaps appear 
in front of these particles, and the faster particles form platoons behind them, 
a phenomenon familiar from vehicular traffic on country roads |TR O, 14" |. The 



system phase separates into regions of density p~ = (the gaps) and regions of 
density p+ = p c (the platoons). 

It is evident from (§) that the condition for phase separation translates into a 
condition on the behavior of the disorder distribution f(p) near p = c. Introducing 
an exponent n through 

f{p)~(p-c) n , p-+c, (5) 

phase separation occurs iff n > 0. At the critical point p = p c the disorder 
averaged headway distribution has a power law tail ~ M^( n + 2 ) [pQ[. Evans 0, |2| 



has emphasized the close analogy to Bose-Einstein condensation, where f(p) plays 
the role of a density of states, and the slowest particle in the system corresponds 
to the quantum mechanical ground state. 

2.2 Coarsening behavior 

No exact results pertaining to the dynamics of phase separation are available, apart 
from the observation [^3] that the existence of a well-defined hydrodynamic limit 



implies that inhomogeneities are restricted to scales smaller than t, and therefore 

lim = 0. (6) 



t- 



Considerable evidence has however accumulated in favor of the idea |20| that the 



coarsening behavior for particlewise disorder can be described in terms of a simpler, 
deterministic model, in which particles move ballistically on the real line with fixed 
random speeds and coalesce upon overtaking. Such a model was first introduced 
by Newell fll3| , and later a detailed kinetic theory was worked out by Ben-Nairn, 



Krapivsky and Redner Q14 



Within the deterministic model, the dynamic exponent z can be determined 
through a simple extremal statistics argument. The key idea is that the particles 
heading the platoons at time t are those with the smallest speeds among of the 
order of £(£) particles. Elementary probability theory suffices to show that, for a 
probability density behaving as ([|), these extremal speeds cluster in an interval 
of size £ _1 /( n+1 ) above the minimal speed c. Therefore the speed difference Av 
between two platoons is of the order ^~ 1 /( ra+1 ) ) and the faster platoon will merge 
with the slower one on a time scale t ~ £/Aw ~ ^( n+2 )/( ra + 1 ). Inverting this relation 
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one obtains the coarsening law ([[]) with 



n + 2 
n + 1' 



(7) 



Numerical results supporting (|7]) have been reported for models with parallel up- 
date p6 |, in simulations of jam dissolution p3 and in a simulation study of a 



system with open boundaries [27 . 



3 Sitewise disorder 
3.1 Disorder types 

We distinguish three cases which will turn out to represent different classes of coars- 
ening behavior. For type I disorder the dynamics is totally asymmetric, q(x) = 0, 
and the forward rates p(x) are independent random variables in an interval [c, 1], 
with a minimal rate c > 0. The simplest (and typical) example is that of binary 
rates, with probability density 

/(p) = 0<J(p-c) + (l-^(p-l) (8) 

where G (0, 1) denotes the fraction of slow sites. Type II disorder is similar to 
type I except that the support of the probability density f(p) extends all the way to 
p = 0, i.e. the minimal rate c = 0. For type II disorder nontrivial dynamics occurs 
only for continuous f(p). As in the models with particlewise disorder (eq.(|5|)), the 
important feature of f(p) is the behavior near p = 0, which can be characterized 
by an exponent n through the relation 

f(p)~p n , p^0. (9) 

Finally, for type III disorder not only the strength, but also the direction of 
the bias is spatially random. A majority of sites has a bias to the right, say, 
with p(x) > q(x), while a minority has q(x) > p(x). If the one-dimensional 
lattice is viewed as a transport path in a higher- dimensional disordered structure, 
such as a percolation cluster, the stretches of minority sites can be interpreted 
as "backbends" where the path turns back against the direction of the driving 
field P8fl . Compared to the strong disorder effects induced by the backbends, the 
randomness in the strength of the bias is irrelevant. Therefore a representative 
example of type III disorder is a model where the strength of the bias is constant, 
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and only its direction varies. This corresponds to setting q{x) = 1 — p(x) and 
choosing the p(x) from a binary distribution which is symmetric around p = 1/2, 



f(p) = (1 - 4>)8{P ~b) + <t>5{p - (i - &))■ (io) 

Here b G (1/2, 1) denotes the strength of the bias and G (0, 1/2) the fraction of 
minority sites. 

It is easy to see that for type II and III disorder the stationary particle cur- 
rent vanishes in the infinite system limit, due to the existence of arbitrarily large 
stretches of arbitrarily small jump rates (for type II) or arbitrarily long backbends 
(for type III). As a consequence phase separation occurs at any density p G (0, 1), 
i.e. p~ = and p+ = 1. For type I disorder the existence of a nontrivial current 
function J(p) > describing the large scale dynamics of density profiles has been 
rigorously established, and it has been shown that J(p) is convex in the sense that 
J"{p) < P9 |. However, in contrast to the models with particlewise disorder the 



stationary state is not known, and therefore an explicit computation of J(p) is not 
possible. In the next section some bounds on J(p) will be derived and used to 
bound the critical densities for type I disorder. The coarsening dynamics for all 



three cases will be addressed in Section 3.3. 



3.2 Bounds on the critical density for type I disorder 

We first collect some obvious properties of J(p). Due to particle- hole symmetry 
we have J(p) = J(l — p). The current is bounded from below by the current 
cp( 1 — p) of a pure system with all rates equal to the minimal rate c, and from 
above by the current p(l — p) of the system with all rates equal to unity. A more 
precise upper bound is obtained by observing that in the infinite system there are 
arbitrarily large stretches with rates arbitrarily close to c. The maximum current 
that can be driven through such a stretch is c/4, the maximum value of cp(l — p). 
We conclude that 



cp(l -p)< J(p) < min[c/4,p(l - p)}. 



Numerical simulations of site-disordered exclusion models [E4l, [L6L 25] and re 



lated growth models |3y, |31| indicate that the upper bound c/4 is attained in a 
finite density interval around p = 1/2, which coincides with the phase separation 
interval [p^T,p^]; by particle-hole symmetry p~ = 1 — p+ = p c . In the following 
our strategy will be to derive optimal lower and upper bounds J<(p) and J>(p) 
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on the stationary current, which are then translated into lower and upper bounds 
pf, on p c through the relation 

J>{ P f) = J<{p > c )=c/A. (12) 

The lower current bound in (JTT|) does not give rise to any nontrivial density bound, 
while the upper bound p(l — p) yields 



p c > (1 - v / T 3 c)/2. (13) 

For the case of binary disorder (eq. (|j) ) an improved lower bound on the current 
was derived by Tripathy and Barma |24| by considering a finite ring of L sites, 
N = pL particles and N s = cf)L slow sites. They start from the observation 
that the maximum current that can be driven through a stretch of slow sites is 
a decreasing function of the length of the stretch (we will return to this point 
below in Section pO|). It is therefore plausible (though not rigorously established) 
that for given L, N and N s the stationary current will be minimal in the fully 
segregated limit where all slow sites form a single large stretch. For L — > oo the 
fully segregated system can be treated as two connected homogeneous systems 
with different densities, which are fixed through the constraints of equal currents 
and total particle number. This yields the upper density bound 



Pc < (1 - (1 - <p)VT^c)/2. (14) 



In the dilute limit — > the bounds dl3|) and fli~4T ) coincide, and give p c = (1 



V 1 — c) /2 exactly. It should however be noted that this limit does not correspond 
to the case of a single defect site, since the maximal current that can be driven 
through a single defect is larger than c/4 p2| (see also Section [3.3|) . 



The lower bound ( Jl3| ) can be improved by comparing the disordered exclusion 
model to a zero range process (ZRP) with the same set of jump rates {p(x)}. 
In the ZRP an arbitrary number of particles is allowed on any site |TJ, [TJ|, and 
therefore any attempted jump succeeds. As a consequence the stationary state 
of the ZRP is a product measure, with the occupation numbers at different sites 
being independent, for any choice of jump rates depending on the position x and 



on the number of particles at the site [19, 33 . Here we consider the case where 



the rate at which a particle is transferred from site x to x + 1 is equal to p(x) 
independent of the number of particles at x, provided the latter is not zero. It is 
then obvious (and can be proved through waiting time considerations) that the 
particle current Jzrp(p) of the ZRP provides an upper bound to the current J(p) 
of the ASEP. 
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In fact the disordered ZRP is equivalent to the ASEP with particlewise disorder, 



with the ZRP occupation numbers representing the headways in the ASEP [19|, [15 
The ZRP current is equal to the particle speed v of the ASEP, which is given 
by (^]) for any disorder distribution f(p). The ZRP density is equal to the mean 
headway of the ASEP, and is therefore related to the ASEP density through pzrp = 
1/Pasep — 1- Evaluating the integral in (fD for the binary distribution @ yields 

Pzrp = Jzbp ; H --. ; , (15) 

\ c — -'ZRP J- — J ZRP J 

and setting J> = Jzrp in ( |P2D we obtain the density bound 

C(1~0) 

Pc > 3 + -j^-, (16) 

which improves (|i3|) for small c. In particular, for c — > we have p c > 0/3, which 
proves, remarkably, that the homogeneous phase p < p c persists even when the 
slow sites become complete blockages. In Figure 3 the bounds ([13]), fll4|) and fll6|) 
are compared to numerical data. 



3.3 Coarsening behavior 

At least for type I disorder the existence of a hydrodynamic limit PD[| implies 
that the relation @ carries over to the sitewise case. To obtain a finer estimate 
of the coarsening scale ^(t) we rely on extremal statistics arguments similar to 
those used in Section |2.2| . A schematic phase separated density profile is shown 
in Figure 4. Two "antischocks" at positions X\ and X2, where the density jumps 
from p+ = 1 — p c to pi = p c , mark bottleneck regions of particularly slow rates, 
which support maximum currents j\ and j%. If the bottleneck in the downstream 
direction is slightly more restrictive, in the sense that Aj = j\ — ji > 0, then the 
low density region between the bottlenecks will slowly fill in and disappear at a 
time 

t « (1 - 2p c )£/Aj. (17) 

If the statistics of extremal bottlenecks is known, the typical current difference Aj 
can be estimated as a function of £ and (0) yields a prediction for the coarsening 
law In the following this will be carried out for the different disorder types. 



3.3.1 Type I disorder 

Consider first the conceptually simplest case of the binary disorder distribution (||). 
We expect long stretches of slow sites to constitute the most restrictive bottlenecks. 
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For a quantitative analysis we would require the maximum current j m ax(c, £) which 
can be driven through a stretch of i slow sites with jump rates c embedded in an 
infinite system of sites with jump rates 1. Already for I = 1 the computation of 
jmax(c, i) is a difficult unsolved problem f32]] . However for large t we can make 
progress by replacing the stretch by a finite system of I sites with uniform jump 
rates c and periodic or open boundary conditions, for which the maximum current 
is known |9|, |10[| . For both kinds of boundary conditions the current approaches 
the I — > oo limit c/4 from above, with a leading correction |34j proportional to 
1/1. Thus we expect, for large £, 

JmaxM) « (c/4)(l + a/£) + 0(1/ f), (18) 

where a is a positive constant of order unity. 

Since the probability distribution of the lengths of slow stretches is 

P(£) = (1 - fitf, (19) 

the longest stretch in a region of size £ is of the order of 



ln(l/0)' 



(20) 



Note that £ max <C £, which is consistent with the assumption of well-localized 
bottlenecks inherent in Figure 4. Using (|T8| ) we see that the currents supported by 



the longest stretches exceed c/4 by an amount of the order of c/£ max , and therefore 

A; « (21) 

In 4 

Inserting this into flTT| ) the leading order coarsening law is obtained as 

with a characteristic time scale to ~ (1 ~ 2p c )/cln(l/0). This argument was 
formulated earlier in the context of phase-disordered growth models, where also 
numerical evidence in favor of the coarsening law (|22|) was presented [^, |3l| . 



For continuous disorder distributions f(p) the identification of the relevant 
bottlenecks is a little more subtle. Consider a region of size t where all rates 
satisfy p(x) < c + e. The maximum current through such a region can then be 
estimated as 

Jmax » ^ (l + jj « c/4 + e/4 + ca/4£ = c/4 + j. (23) 
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If f(p) behaves as in (||) for p — > c the probability of the region is of the order of 
e (n+i)i_ rjj ie probability distribution of j can then be written as 



P(j) ~ Jdijde e (n+1) %' - (e + ca/l)/A) ~ 

J deexp[-(n + l)caln(l/e)/(4j-e)]. (24) 

Evaluating the last integral at the saddle point yields 

P(j) ~ exp[-(n + l)caln 2 (l/j)/j']. (25) 

The current scale Aj of the most restrictive among £ bottlenecks is obtained by 
setting P(Aj) ~ l/£, which gives 

A, ~ ^ (26) 

and the corresponding coarsening law reads, to leading order in t, 

t\\n(\nt)} 2 . . 

which, for most purposes, is indistinguishable from (p2[) . 
3.3.2 Type II and III disorder 

For type II disorder with a continuous probability distribution f(p), characterized 
by (H), the expression (p3|) for the maximum current supported by a slow stretch 
of length £ applies with c = 0. The distribution of j max then becomes 

P(jmax) -Jdijde e(" +1 )% max - e(l + a/i)/A) ~ 

d£exp[-(n + l)£ln((l + a/£)/ij max )}. (28) 

Now the maximum of the exponent evidently occurs at I = 1, i.e. the dominant 
bottlenecks are individual slow sites. The distribution of the currents supported 
by the bottlenecks is then simply given by the jump rate distribution f(p) itself, 
and the situation reduces to that analyzed in the case of particlewise disorder, 
Section |2.2| . In particular, the coarsening exponent z for type II sitewise disorder 
is also given by (0). 
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For type III disorder with distribution (|H]) the dominant bottlenecks are long 
backbends, i.e. stretches of minority sites at which the local bias is directed against 
the mean flow direction. The maximum current that can be driven through a 
backbend of length I is exponentially small in £, and is given by |24], |28 



JmaxW ~ exp[-(l/2)nn(6/(l - 6))]. (29) 

Combining this with the probability distribution (|19"D of backbend lengths it follows 
that j max is distributed according to a power law, 

^(Jmax) ~ (jmax) 29 "" 1 (30) 

where 

wm ■ (3 ) 

Since the largest backbend in a region of size £ is of length £ ~ ln£ <C £, we can 
employ a coarse grained picture in which the backbends are shrunk to individual 
sites with a jump rate distribution given by (|3~0D , thus effectively reducing the 
problem to type II disorder with the exponent n in (|j) given by n = 2/8 — 1. The 
coarsening exponent for the disorder distribution ( [T0|) is then obtained from (]?]) as 

z = 1+9/2. (32) 
3.4 Relation to directed polymers 

Using the waiting time approach |35|] the site disordered ASEP can be mapped to 
a zero temperature directed polymer (DP) with point and columnar disorder ||36|| . 
In that context the coarsening law £(t) describes the disorder-induced transverse 
wandering of the polymer, which can be estimated using variable range hopping 



arguments [35 and the analogy to Lifshitz tails for one-dimensional disordered 



Schrodinger operators ]37 



To see that the results derived for the DP are consistent with those obtained 



above, it is important the recall [|35| that the waiting time mapping transforms 



the time t of the ASEP into the energy of DP. For type I disorder the transverse 



wandering 5x of the DP was found to increase with its length L as [|36 



5x~L/(lnL) 2 , (33) 

while the ground state energy behaves as E ~ L/ In L to leading order. Combining 
the two results and identifying E ~ t the coarsening law (E2I) follows. 



11 



For type II disorder the power law @ of the probability distribution at small 
p translates into a power law tail 



P(r) 



~ r 



-(2+ra) 



OO 



(34) 



in the distribution of waiting times or energies r = I /p. Directed polymers in 
the presence of columnar disorder with a power law distribution were considered 
in Ref.[37 , where it was shown that the wandering is typically ballistic, 5x ~ L, 
while the ground state energy scales with length as 



E ~ L {n+2),{n+l) 



(35) 



in agreement with ([/]). It is worthwhile to point out that in the DP context the 
scaling laws (|33l) and (35) were also confirmed numerically [36, 371. 



4 Summary and open questions 

In this paper I have described some recent progress in our understanding of dis- 
order effects in asymmetric simple exclusion models. A common feature of both 
particlewise and sitewise disordered systems is the appearance of phase separation 
in an interval of densities, which is macro scopically characterized by a linear por- 
tion in the current-density relation J(p); in the particlewise case J(p) = cp for 
p < p c , while in the sitewise case J(p) = c/4 for p c < p < 1 — p c . An interest- 
ing open question concerns the connection between phase separation and linearity 
of J(p), which is reminiscent of the role that the convexity of thermodynamical 
potentials plays for the stability of equilibrium systems. While it is obvious that 
phase separation implies a linear segment in J(p), the converse statement has, to 
my knowledge, not been established. To prove that it is false, it would be sufficient 
to find a (noisy!) exclusion type model with a homogeneous stationary state and 
a linear current-density relation (deterministic systems with linear J(p) are well 
known |7|,|3§). 

The dynamics of phase separation has been explored in the framework of scaling 
arguments, which can be formulated in a similar way both for particlewise and 
sitewise disorder. In the particlewise case the relevant bottlenecks which determine 
the positions of domain boundaries are always individual slow particles, while in 
the sitewise case with type I and III disorder the bottlenecks are formed collectively 
by many defects. For type I disorder this implies a certain universality of the 
coarsening law, in the sense that the exponent z in ([1]) is z — 1 independent of the 
underlying disorder distribution; the additional logarithmic corrections in (|22| , |2TD 
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ensure the consistency with the rigorous result (||). This is somewhat analogous 
to the case of finite temperature directed polymers with columnar defects, where 
universal scaling laws arise from the thermal averaging over large spatial regions 



A numerical confirmation of the predictions for the coarsening dynamics in the 
case of sitewise disorder would be most welcome. For type II disorder this should 
be relatively straightforward, however in the cases of type I and III disorder the 
behavior is dominated by exponentially rare regions, which may make it hard to 
reach asymptopia. 
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Figure 1: Illustration of the asymmetric simple exclusion process. 




Figure 2: Space-time plot of trajectories of the ASEP with particlewise disorder. 
The figure shows 256 particles on a ring of 1024 sites. The initial distribution of 
particles was random with density p = 1/4 < p c = 0.4. Courtesy of M. Gerwinski. 
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Figure 3: Bounds on the critical density for sitewise disorder with the binary 
distribution = 1/2. The dashed line is the upper bound ([TID , the two full 
curves show the lower bounds (|13|) and (|T6|), and the full squares are simulation 
data obtained by G. Tripathy. 
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Figure 4: Schematic of two bottlenecks at positions X\ and x 2 . 



